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ABSTRACT 

We have investigated the molecular environment of the semicircular composite supernova remnant 
(SNR) 3C 396 and performed a Chandra spatially resolved thermal X-ray spectroscopic study of 
this young SNR. With our CO millimeter observations, we find that the molecular clouds (MCs) 
at Vlsr ~ 84 km s _1 can better explain the multiwavelength properties of the remnant than the 
Vlsr= 67-72 km s _1 MCs that are suggested by Lee et al. (2009). At around 84kms _1 , the western 
boundary of the SNR is perfectly confined by the western molecular wall. The CO emission fades 
out from west to east, indicating that the eastern region is of low gas density. In particular, an 
intruding finger /pillar-like MC, which may be shocked at the tip, can well explain the X-ray and radio 
enhancement in the southwest and some infrared filaments there. The SNR-MC interaction is also 
favored by the relatively elevated 12 CO J=2-l/J=l-0 line ratios in the southwestern "pillar tip" and 
the molecular patch on the northwestern boundary. The redshifted 12 CO (J=l-0 and J=2-l) wings 
(86-90 km s _1 ) of an eastern 81 krns" 1 molecular patch may be the kinematic evidence for shock-MC 
interaction. We suggest that the 69 km s _1 MCs are in the foreground based on HI self- absorption 
while the 84 km s _1 MCs at a distance of 6.2 kpc (the tangent point) are in physical contact with 
SNR 3C 396. The X-ray spectral analysis suggests an SNR age of ~ 3 kyr. The metal enrichment of 
the X-ray emitting gas in the north and south implies a 13-15M B1-B2 progenitor star. 

Subject headings: ISM: individual (3C 396, G39.2— 0.3) - ISM: molecules - supernova remnants 



1. INTRODUCTION 

The progenitors of core-collapse supernovae are mas- 
sive stars and they most likely formed in giant molec- 
ular clouds (MCs). Because of their short lifetime, one 
may expect that the supernova remnants (SNRs) of core- 
collapse supernovae would be near their natal environ- 
ment. The association with MCs often makes SNRs dis- 
play irregular morphologies in multiwavelengths, which 
indicates complicated shock interaction with the com- 
plex, inhomogeneous ambient interstellar medium (ISM). 

3C 396 (G39.2— 0.3) is a semicircular young compos- 
ite SNR, exhibiting both the central X-ray bright pulsar 
wind nebula (PWN) and western incomplete half-shell 
structures in multiwavelengths. In the high dynamic- 
range VLA images of the remnant, a remarkable "tail" 
extending from the very faint eastern shell of the SNR 
is revealed and is interpreted as an eastern "blow-out" 
structure (Patnaik et al. 1990). From VLA 6 and 20 cm 
data, a large scale brightness asymmetry is found in the 
east-west direction across SNR 3C 396, the western half 
of the remnant having roughly twice the flux density of 
the east (Anderson & Rudnick 1993). The Einstein X- 
ray emission is bright in the limb and the center of the 
remnant (Becker & Helfand 1987). Based on ASCA X- 
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ray spectral analysis, a nonthermal component in the 
center of 3C 396 is attributed to a synchrotron PWN 
(Hurrus & Slane 1999). Recently, with about 100 ks 
high spatial resolution Chandra observation, a nonther- 
mal nebula containing a point-like source is also found 
in the center of the remnant and is interpreted as a 
synchrotron PWN surrounding a yet undetected pulsar 
(Olbert et al. 2003). The possibility of the excess mi- 
crowave emission arising from a spinning dust compo- 
nent toward 3C 396 is investigated at 33 GHz with the 
Very Small Array (VSA) telescope (Scaife et al. 2007). 
Based on the detection of H2 molecules and the sharp 
southwestern boundary in the radio emission, the rem- 
nant is suggested to be encountering a MC in this region 
(Lee et al. 2009). The remnant is suggested at a dis- 
tance of about 8.5 kpc because a cavity-like structure of 
the Vlsr ~ 69 km s _1 molecular gas seems to surround 
the remnant. The Spitzer IRS observation reveals that, 
along the western shell, the inner [FeII]-line filament and 
the outer H2-line one is found to be spatially separated 
(Hewitt et al. 2009), which is similar to the result of 
near-IR 1.64 /im [Fell] and 2.12 /im H2 line observations 
using Wide-field Infrared Camera aboard the Palomar 5 
m Hale telescope (Lee et al. 2009). It is most probably 
because of the various molecular gas density conditions 
of the multi-phase preshock ISM across the western shell 
of 3C 396, where the shock of the remnant is interact- 
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ing with a complicated distribution of MCs (Hewitt et 
al. 2009). A detailed investigation of the molecular gas 
environment and the hot gas interaction with it are thus 
necessary, which may be helpful to understand the mul- 
tiwavelength physical properties and the dynamic evolu- 
tion of this irregularly shaped SNR. 

In this paper, we have presented our millimeter CO line 
observations toward SNR 3C 396 and a spatially resolved 
Chandra thermal X-ray analysis. The observations and 
the data reduction are described in Section 2. The main 
observational results and calculations are presented in 
Section 3. The discussion and conclusions are given in 
Section 4 and Section 5, respectively. 

2. OBSERVATIONS AND DATA REDUCTION 

The observation was first made in the 12 CO (J=l- 
0) line (at 115.271 GHz) in 2006 April using the 6 m 
millimeter-wavelength telescope of the Seoul Radio As- 
tronomy Observatory (SRAO) with single-side band fil- 
ter and the 1024-channel autocorrelator with 50 MHz 
bandwidth. The half-power beam size at 115 GHz is 2' 
and the typical rms noise level was about 0.2 K at the 
0.5 km s _1 velocity resolution. The radiation temper- 
ature is determined by Tr = Tx/{fh x to), where Ta 
is the antenna temperature, /b the beam filling factor 
(assuming /b ~ 1), and 7? m b the main beam efficiency 
(~ 75%). We mapped the 14' x 15' area covering 3C 396 
centered at (19 h 04 m 06 s , 05°27'00") with grid spacing of 

The follow-up observations were made simultane- 
ously in the 12 CO (J=l-0) line (at 115.271 GHz), 
13 CO (J=l-0) line (110.201 GHz), and C 18 (J=l-0) 
line (109.782 GHz) during 2008 November-December us- 
ing the 13.7 m millimeter- wavelength telescope of the 
Purple Mountain Observatory at Delingha (hereafter 
PMOD). The spectrometer has 1024 channels, with a 
total bandwidth of 145 MHz (43 MHz) and the veloc- 
ity resolution of 0.5 km s" 1 (0.2 km s" 1 ) for 12 CO( 13 CO 
and C 18 0). The mean rms noise level of antenna temper- 
ature was less than 0.1 K for those lines. The intensity 
of the C 18 (J =1-0) emission is too weak to examine 
its spatial distribution. The half-power beamwidth of 
the telescope is about 54" and the main beam efficiency 
^mb was 52%. We mapped the 9' x 9'. 5 area centered at 
(19 h 04 m 03 s ,05°26 / 30 // ) with grid spacing of 30". 

We also made the 12 CO (J =2-1) line observation to- 
ward partial regions of the SNR with grid spacing of 
30" in 2010 February using the telescope of the SRAO 
to complement the data. The half-power beam size at 
230 GHz is 48", the main beam efficiency was ~ 57%, 
and the typical rms noise level was about 0.2 K at 
the 0.13 km s _1 velocity resolution. The resolution of 
the SRAO 12 CO (J=2-l) observation and the PMOD 
12 CO (J =1-0) observation are very similar, which makes 
it convenient to directly compare the two observations. 

All the CO observation data were reduced using the 
GILDAS/CLASS package developed by IRAM0 

The Chandra X-ray and Spitzer 24 /im mid-IR ob- 
servations are also used for analysis. We revisited the 
Chandra ACIS observational data of SNR 3C 396 (Ob- 
sIDs: 1988; PI: S. P. Reynolds) with total exposure time 

2 http://www.iram.fr/IRAMFR/GILDAS 



of about 100 ks. We reprocessed the event files (from 
Level 1 to Level 2) using the CIA03.4 data processing 
software to remove pixel randomization and to correct 
for CCD charge-transfer inefficiencies. The overall light 
curve was examined for possible contamination from a 
time-variable background. The reduced data, with a net 
exposure of 93.7 ks, were used for our subsequent anal- 
ysis. The mid-IR 24 /im observation used here was car- 
ried out as the 24 and 70 Micron Survey of the Inner 
Galactic Disk Program (PID: 20597; PI: S. Carey) with 
the Multiband Imaging Photometer (MIPS) (Rieke et 
al. 2004). The Post Basic Calibrated Data (PBCD) of 
the 24 /im mid-IR were obtained directly from Spitzer 
archive. The VLA 1.4 GHz radio continuum emission 
(L band) data were adopted from Anderson & Rudnick 
(1993). The Hl-line data and the large scale 13 CO line 
data were obtained from the archival VLA Galactic Plane 
Survey (VGPS: Stil et al. 2006) and the Galactic Ring 
Survey (GRS: Jackson et al. 2006), respectively. 

3. RESULTS 

3.1. Molecular Gas Environment of SNR 3C 396 

The large field of view (14' x 15') 12 CO channel maps 
are made with the SRAO observation and shown in 
Figures [1] and [2j overlaid with radio (red thick) and 
X-ray (cyan thin) contours. The velocity components 
around 69 km s _1 (in 64-72 km s _1 interval) and around 
84 km s _1 (in 80-88 km s _1 interval) both seem to have 
positional correlations with SNR 3C 396. The 69 km s _1 
MCs are mainly situated in the north of the remnant 
[with the 64-67 km s _1 molecular gas in the northeast 
and the 68-71 km s _1 gas in the northwest (NW)]. The 
SNR seems to be projected roughly within a cavity-like 
structure of the 12 CO 69 km s _1 component (Figured]). 
However, the morphology of the 69 km s _1 component 
appears to not well match the radio boundary of the rem- 
nant; especially, the 68-71 km s _1 gas appears to over- 
lap the northwestern corner of the remnant. On the other 
hand, at around 84 km s _1 (Figure [2]), a thick wall struc- 
ture of molecular gas is seen in the west; the eastern face 
of the wall at the 85-87 km s _1 interval perfectly follows 
the western radio boundary of the remnant. The radio 
contours of the remnant seem to coincide with a cavity 
at 85-87 km s _1 . The CO emission around 84 km s _1 
and the remnant's radio emission both seem to fade out 
gradually from west to east. 

These characteristics are seen more clearly in the close- 
up image (Figure [3j left panel) made from the PMOD 
observation with a higher resolution. As also shown in 
Figure [3] (left panel), the X-ray/mid-IR/radio emission in 
the west of the remnant appears to be well confined by 
the western molecular wall around 84 km s _1 . Moreover, 
a prominent finger/pillar-like MC is seen in the southwest 
(SW), with one end intruding inside the SNR border. 
The radio and IR brightness peaks of the remnant are 
located at the short southwestern section of shell that 
"cuts" the molecular pillar. Also in the southwestern 
edge, an X-ray enhancement is coincident with the radio 
peak; in Section 13^21 below, we will show that this X-ray 
enhancement is exceptionally soft and of normal metal 
abundance. 

As an attempt to remove the confusion by large-scale 
density gradient, we used the unsharp masking method 
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as described in Landecker et al. (1999), in which we sub- 
tracted the smoothed 12 CO intensity map from the orig- 
inal 12 CO map. Similar images can also be produced by 
subtracting the smoothed and scaled 13 CO intensity im- 
age from original 12 CO image. In the "cleaned" image 
of the 69 km s -1 MCs, the cavity-like structure is seen 
again (Figure HJ the right panel in the first line), but 
it clearly does not well match the radio boundary of the 
remnant. For the 84 km s -1 MCs, on the other hand, the 
southwestern finger/pillar-like intrusion and the western 
molecular wall can be seen clearer in the "cleaned" im- 
age (Figure IH the right panel in the second line) than in 
the intensity maps. The western boundary of the SNR 
appears to perfectly follow the left side of a molecular 
wall again (Figure |4j the right panel in the third line). 
A patch of strong 12 CO emission around 84 km s -1 is 
seen in the NW, where the radio emission of the SNR 
is slightly enhanced; this patch can also be discerned in 
the channel map, Figure [2j at 84 kms -1 . In addition, 
there is another patch of 12 CO emission in 87-90 km s -1 
interval in the eastern "blow-out" region of the remnant 
(Figure HI the right panel in the fourth line); it can also 
be discerned in Figure [2] at 87-88 km s -1 . 

The integrated C O-line intensity ratio between the dif- 
ferent transitions is suggested as a probe of the SNR- 
MC interaction (Seta et al. 1998; Jiang et al. 2010). 
The 12 CO J= 2-1/ J= 1-0 ratio maps towards 3C396 
at 66.5-68, 68-70.5, 70.5-72, 80-83, 83-85, and 85- 
88 kms -1 (which covers the centers and wings of the 
69 and 84 kms -1 lines, respectively) are presented in 
Figure [5] For the ~ 84 kms -1 component, the ratios 
are clearly elevated to as high as > 0.9 in 80-83 kms -1 
interval (compared with ~ 0.5-0.7 in average) in the 
southwestern molecular intrusion region and the north- 
western molecular patch (corresponding to the two radio 
brightness peaks, respectively) described above. For the 
~ 69 kms -1 (66.5-72 kms -1 ) component, although the 
ratios are weakly elevated in some boundary regions, the 
ratios are all below 0.9. Therefore the elevated 12 CO J= 
2-1/ J= 1-0 ratios seem to better imply an interaction 
of the SNR with the ~ 84 kms -1 MCs than with the 
~ 69 kms -1 ones. 

We have carefully examined the CO spectra in the field 
of view and presented the 12 CO( J = 1-0), 12 CO( J = 2- 
1), and 13 CO(J = 1-0) spectra extracted from two re- 
gions (in the east and SW, as marked in the right panel 
in the fourth line of Figure 2J in Figure [6] Around 
69 kms -1 , the 13 CO line profiles are very similar to the 
12 CO profiles for both the "East" and "SW" regions. 
Both the 12 CO and 13 CO spectra of the "SW" region 
have peaks at 84 km s -1 , with broad blue (left) wings 
in the line profiles. Yet, it is difficult to judge whether 
the blue wings of the 84 km s -1 12 CO lines are dynami- 
cally broadened or not, due to the significant blue wing 
of the 13 CO line and especially a distinct 13 CO peak 
at 81 kms -1 . 13 CO emission is usually optically thin 
and indicative of quiescent, intrinsically high-column- 
density MCs. For the "East" region, both the 12 CO 
and 13 CO spectra have peaks at 81 km s -1 ; notably, 
there is an additional broadening in the red (right) wing 
(86-90 km s -1 ) of the 81 km s -1 12 CO line, as com- 
pared with little 13 CO emission in the same velocity in- 



terval. It most likely indicates that the 12 CO emission in 
this interval results from the perturbation suffered by the 
81 km s -1 molecular gas in the east of the remnant. Al- 
though there is a velocity discrepancy between 81 kms -1 
for the eastern region and 84 kms -1 for the western re- 
gion, the undisturbed clouds in the SNR field may have 
a wide velocity distribution, e.g., the LSR velocity of 
MCs with which SNR IC 443 interacts is in the range 
of —7 to —2 kms -1 (van Dishoeck et al. 1993; Hewitt 
et al. 2006). Thus the 81 kms -1 cloud in the east is 
very likely, together with the 84 kms -1 clouds, among 
the MCs shocked by SNR 3C 396. 

The above analysis show that both of the 69 and 
84 km s -1 molecular components appear to have spa- 
tial correlation with SNR 3C 396, but the latter seems 
to have much better evidence than the former. We will 
further show that the 69 km s -1 clouds are in the fore- 
ground of SNR 3C 396, using the VGPS radio and GRS 
13 CO data. 

The 1420 MHz continuum image (Figure [7J upper left 
panel) displays semicircularly brightened structure in the 
side roughly close to the Galactic plane and faint emis- 
sion in another side. Following the method used in Tian 
et al. (2007), we produced the HI spectrum (Figure [7J 
upper right panel) of the radio peak region of the rem- 
nant. The maximum radial velocity of absorption fea- 
tures in the direction towards 3C 396 is up to the tangent 
point LSR velocity (i.e., ~ 84 km s -1 ), and virtually no 
significant absorption feature appears at negative veloc- 
ities. It indicates that the distance of SNR 3C 396 is 
between 6.2 kpc (the tangent point) and 12.5 kpc (the 
Solar circle). Here we have used the rotation curve in 
Clemens (1985) with R = 8.0 kpc (Reid 1993) and 
Vq = 220 km s -1 and adopted the spiral model in Taylor 
& Cordes (1993; Figure 1 therein). 

The HI self-absorption (SA) measurement can provide 
another clue to the MCs' distances. A single radial ve- 
locity indicates two distances (the near and the far kine- 
matic distances) in the inner Galaxy along the line-of- 
sight (LOS). MCs in the near distance often display HI 
SA features and many work have been carried out based 
on this property (e.g., Liszt et al. 1981; Jackson et al. 
2002; Goldsmith & Li 2005; Busfield et al. 2006; Ander- 
son & Bania 2009; Roman-Duval et al. 2009). The HI 
SA property also seems effective in resolving kinematic 
distance ambiguity toward SNR-MC associations, such 
as 3C 397 (Jiang et al. 2010) and CTB 109 (Tian et al. 
2010). Figure [7] (lower panels) shows the VGPS HI and 
GRS 13 CO spectra of two regions (as defined in the up- 
per left panel) near the remnant. The HI SA feature is 
clearly seen at 69 km s -1 associated with the 13 CO peak, 
which indicates that the 69 km s -1 molecular component 
is most probably located at the near distance along the 
LOS, i.e., 4.3 kpc. This component seems to be a part of 
the 69.65 km s -1 MC complex GRSMC G039.34-00.26 
at the near side (Roman-Duval et al. 2009), which is ap- 
parently outside the distance range 6.2-12.5 kpc inferred 
for SNR 3C 396. 

Therefore we suggest that SNR 3C 396 is more prob- 
ably associated with the MCs at ~ 84 km s -1 . The 
derived column density and gas mass of some molecu- 
lar features (the northwestern patch, southwestern pillar 
and the tip, and the western wall (Figure |HJ)) are pre- 
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sented in Table [TJ In the derivation of these parameters, 
two methods are used and similar results are yielded. In 
the first method, the H2 column density is calculated 
by adopting the mean CO-to-H2 mass conversion factor 
7V(H 2 )/W( 12 CO)^ 1.8 x 10 20 cm^K^km^s (Dame 
et al. 2001). In the second method, on the assump- 
tion of local thermodynamic equilibrium and the exci- 
tation temperature ~ 10 K, the 13 CO column density is 
converted to the H2 column density using the relation 
7V(H 2 )=7 x 10 5 7V( 13 CO) (Frerking et al. 1982). Assum- 
ing a cuboidal volume for the wall region (Figured]) with 
size 4'. 5 x 4'. 5 x 9'. 5, the mean molecular density of the 
western wall is inferred as ~ 1 70d^ H 2 cm . 

3.2. X-ray Properties of SC 396 

Next, we investigate the positional relation of 
SNR 3C 396 X-ray emission with the CO intensity distri- 
bution in detail. In general, we extracted Chandra, X-ray 
spectra from five X-ray bright regions (as defined in Fig- 
ure [9]) using the CIA03.4 script specextract. The back- 
ground spectrum was extracted from four nearby regions 
(dashed circles in Figure [9]). The five on-source spec- 
tra were adaptively regrouped to achieve a background- 
subtracted signal-to- noise ratio (S/N) of 5 per bin. Then 
we used a nonequilibrium ionization (NEI) model to fit 
the spectra of the remnant. The XSPEC spectral fit- 
ting package was used throughout. For the foreground 
absorption, the cross sections from Morrison & McCam- 
mon (1983) were used. The parameters of the fitting are 
listed in Table [2] with the 90% confidence ranges. 

The five spectra all contain distinct lines of Si XIII 
(~ 1.85 keV) and S XV (~ 2.45 keV), but Ar XVII 
(~ 3.12 keV) and Ca XIX (~ 3.86 keV) lines are promi- 
nent only in the northern and southern regions (Fig- 
ure [10]). The model of our spectral fit generates a vari- 
able hydrogen column density from 4.0 x 10 22 cm -2 to 
6.5 x 10 22 cm -2 . These values are roughly consistent with 
4.65 x 10 22 cm -2 from the diffuse X-ray spectral fitting 
with ASCA X-ray observation (Harrus & Slane 1999) and 
5.3 x 10 22 cm -2 from the entire central nebula spectral fit- 
ting with Chandra X-ray observation (Olbert et al. 2003). 
The northern and southern regions have relatively high 
temperatures (> 1 keV) and slight overabundance of Si, 
S, and Ca. The temperatures of the other three regions 
are relatively low, around 0.7 keV. The western edge is 
also S enriched and, notably, only the southwestern re- 
gion is found to be of normal abundance. The ionization 
timescales, n e t, in the northern, southern, and western 
regions are less than 5.5 x 10 11 s cm -3 but seems to be 
higher in the southwestern region (> 4 x 10 11 s cm -3 ). 

An estimate of the hot-gas hydrogen-atom density nx 
of these regions is made using the volume emission mea- 
sure, fn e nxV, derived from the normalization of the 
XSPEC model, where n e is the electron density, V the 
volume of the emitting region, and / the volume filling 
factor of the hot-gas (0 < / < 1). Here we have assumed 
n e ~ 1.2nx and oblate spheroids for the defined elliptical 
regions (Figure EJ). The derived physical parameters are 

listed in Table [3] If we adopt nx = 2f~ 1 / 2 d 6 ^ 2 cm -3 
as the mean density of SNR's hot-gas (where 6^.2 = 
d/6. 2 kpc is used for scaling; Section l47Tj) and assume an 
oblate spheroid (2'. 7 x 3'. 9 x 3' .9) for the entire diffuse 
thermal X-ray emitting volume, V = 3.0xl0 58 fd\ 2 cm 3 , 



the mass of the X-ray emitting hot-gas of SNR 3C 396 is 

In the Chandra ACIS three-color image (Figure [3] right 
panel), the X-ray emission is hard in the central PWN 
region and the northern and southern regions (domi- 
nated by the blue color), relatively soft along the western 
boundary (dominated by the green color), and softest in 
the southwestern edge (seen as an orange patch). This is 
consistent with the distribution of hot-gas temperatures 
derived above. The softest patch is coincident with the 
tip of the Vlsr ~ 84kms _1 molecular finger/pillar. 

4. DISCUSSION 

4.1. Association of SNR SC 396 with the 
Vlsr ~84 km 5 _1 MCs and the Distance 

The cavity-like 12 CO structure at Vlsr ~ 69kms _1 
seems consistent with the GRS 13 CO data (Lee et al. 
2009; Figure 7 therein). Based on such a structure at 
^lsr ~ 69kms _1 that seems to roughly surround SNR 
3C 396, Lee et al. (2009) suggest an association of 3C 396 
with the ~ 69kms _1 MCs and thus adopt a distance 
8.5 kpc. However, the HI SA property of the 69 km s _1 
molecular component indicates that it is most likely to be 
in the foreground (~ 4.3 kpc) toward SNR 3C 396, while 
the 21 cm HI absorption feature suggests the SNR bein g 
located in the distance range 6.2-12.5 kpc (Section l3TTj) . 
It thus appears difficult to associate the ~ 69 kins" 1 MCs 
with the remnant. In our 12 CO observation, we indeed 
find that the cavity-like structure of the 69 km s _1 molec- 
ular component does not well match the morphology of 
the remnant (Section 13. ip . 

Besides the 69 km s _1 MCs, however, we have found 
that the western boundary of the SNR is perfectly 
confined by the western molecular wall at Vlsr ~ 
84kms _1 and the multiwavelength (X-ray, mid-IR, and 
radio) properties are consistent with the presence of the 
84kms- 1 MCs (Section GO}. 

The eastern 12 CO patch (seen in Figure [2] at 87- 
88kms _1 ; Figured] the right panel in the fourth line) is 
shown to have a redshifted broadening (87-90 km s _1 ) in 
the 12 CO spectral line profiles (Figure [6] left panel) and 
thus probably represents the molecular gas perturbed by 
the SNR. The molecular line broadening has widely been 
used as a kinematic evidence of interaction between SNR 
shocks and MCs (see Jiang et al. 2010), such as in the 
studies of prototype interacting SNRs IC 443 (DeNoyer 
1979), W28 (Arikawa et al. 1999), HB21 (Byun et al. 
2006), and G347.3-0.5 (Moriguchi et al. 2005). More re- 
cently, broadened line profiles of optically thin 12 CO lines 
were also successively found in our serial work on SNR- 
MC interaction (Kes 69, Zhou et al. 2009; Kes 75, Su et 
al. 2009; and 3C 397, Jiang et al. 2010). 

Although the broad blue wings of the ~ 84kms _1 
12 CO lines of the western region of 3C 396 cannot be ex- 
clusively determined to be dynamically broadened pro- 
files (Section 13. 1( Figure [6] right panel), the relatively 
enhanced 12 CO J= 2-1/ J= 1-0 ratios in such wings (80- 
83kms _1 ) in the northwestern and southwestern bound- 
aries (Figure [5] lower left panel) coincide with the south- 
western molecular intrusion region and the northwest- 
ern molecular patch, and also with two radio brightness 
peaks, respectively, which is in favor of a perturbation 
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of the western molecular gas by the SNR expansion. We 
note that the maximum ratio obtained here is not higher 
than 1, while higher values are observed in other in- 
teracting SNRs (e.g., W44, IC 443, HB21, G349.7+0.2, 
and G18.8+0.3, see Dubner et al. 2004 and references 
therein). Two main reasons may explain why the ra- 
tio value in 80-83 km s _1 interval does not exceed unity: 
1. The filling factor of the disturbed 12 CO (J=2-l) emit- 
ting region may be small in the 48" beam size in 3C 396 
compared to that for 12 CO (J=l-0), and the line ratios 
could thus be diluted. 2. The location of 3C 396 at the 
tangent point makes the environment much complicated 
and may cause unignorable velocity crowding in the in- 
terval of interest, with the line wings difficult to differen- 
tiate. Actually, e.g., there is a distinct cloud represented 
by the 13 CO line peak at 81kms _1 in the southwestern 
molecular intrusion region (see Figure 6, right panel). 
Confusion of the 12 CO emission of such gas will decrease 
the 12 CO J= 2-1/ J= 1-0 ratios in the velocity interval. 

In multiwavelength morphologies, both the X-ray and 
radio emissions are bright in the western half and faint 
in the eastern half (Figures [3] and [9j), which is proba- 
bly because of the large-scale density gradient roughly 
perpendicular to the Galactic plane. The bright side of 
the remnant is close to the Galactic plane (see Figure [7J 
upper left panel, in Galactic coordinates). The bright 
X-ray, IR and radio emissions in the west are well con- 
fined by the thick molecular wall of the ~ 84 km s _1 
component, which may suggest that the SNR expan- 
sion is hampered there. This scenario is strengthened 
by the radio and IR brightness peaks and the soft X- 
ray patch in the southwestern edge, which are coincident 
with the "tip" of the intruding finger /pillar of molecular 
gas. When the SNR shock hits the dense molecular pillar 
or clumps, it will be decelerated and thus the magnetic 
field be compressed and amplified. This can account for 
the prominent radio enhancement and sharp edge there. 
Actually, the magnetic field in the SW of the remnant is 
found to be essentially tangential to the shell, which is 
probably because the random interstellar magnetic field 
has been chiefly swept up into the accumulated gas and 
compressed along with the remnant's shell (Patnaik et al. 
1990). [Parenthetically, the tangential field (Patnaik et 
al. 1990) and the slightly enhanced radio emission in the 
northwestern edge are consistent with the presence of the 
dense molecular clump in the NW (Section f3TTj) .] In X- 
rays, the low electron temperature, normal metal abun- 
dance, and large ionization timescale (as obtained from 
our spectral fitting) of the soft southwestern X-ray patch 
can be naturally interpreted by the interaction with the 
dense molecular pillar. The IR morphology can also be 
accounted for with the scenario of shock interaction with 
western molecular £ as, as discussed in the following. 

In IR, there is a spatial separation between the inner 
near-IR 1.64 /am [Fell] and outer 2.12 jam H2 emission fil- 
aments along the western edge of the remnant (Lee et al. 
2009; Hewitt et al. 2009). Lee et al. (2009) interprets this 
morphology due to the SNR shock overtaking the wind 
of a red supergiant (RSG), giving rise to [Fell] emission 
and the interaction with a dense circumstellar bubble, 
giving rise to H2 emission. Such a RSG wind bubble 5- 
7 pc in radius invokes a 25-35M progenitor. This mass 
assumption is different from our estimate of 13-15M 



derived from the X-ray emitting metal abundances (see 
below in Section 14. 2. 3p and is not in agreement with the 
mass range < 25M calculated for the neutron star pro- 
genitors with metallicity Z < Z (Heger et al. 2003). 
The 13-1 5 Mq progenitor of type-IIP supernova (see Sec- 
tion I4.2.3|) can only blow a RSG wind bubble of radius 
< 1 pc (Chevalier 2005). However, the spatially sepa- 
rated [Fell] and H2 emission can be interpreted in term 
of SNR-MC interaction. As Reach et al. (2005) and He- 
witt et al. (2009) suggest, the H2 emission can arise from 
the shock interaction with dense clumps while the [Fell] 
emission can arise from the same shock encountering the 
less dense material at the clumps' surface or from the 
eject a. 

The bright 24 /im mid-IR emission (Figure [3]) may 
chiefly come from the dust grains and even possibly from 
the ionic and molecular species in the shocked gas (e.g., 
Tappe et al. 2006; Su & Chen 2008) when the western 
molecular wall is hit by the blast wave. A spinning dust 
component is invoked to account for the spectral excess 
at 33 GHz (Scaife et al. 2007), and it seems that the 
needed dust may lie in the MCs and be heated by the 
shocked gas. 

All these evidences strongly point to a conclusion that 
SNR 3C 396 is associated, and in physical contact, with 
the - 84kms _1 MCs. The LSR velocity 84 km s _1 cor- 
responds to a Galactic location near the tangent point, 
at a distance of ~ 6.2 kpc; this dis tanc e is consistent 
with the range [6.2-12.5 kpc (Section l3.ip ] inferred from 
the HI absorption data. We also notice that a distance 
of 6.5 kpc suggested by Hewitt et al. (2009) based on the 
hydrogen column is very similar to our result. Hereafter, 
we write the distance to 3C 396 as d — 6.2^6.2 kpc for 
scaling. 

It was queried by Patnaik et al. (1990) whether the 
HII complex NRAO 591 (or U39.25-0.07), which looks 
projectively near 3C396, is associated with the SNR. 
Since the complex is at Vlsr= 23 km s _1 with the far 
side distance of 12.1 kpc (Anderson & Bania 2009), it 
is not associated with the remnant. Instead, it is very 
likely to be associated with the 22.46 km s" 1 GRSMC 
G039.24-00.06 at the far side (Roman-Duval et al. 2009). 

4.2. SNR physics of 3C 396 
4.2.1. Global evolution 

According the above multiwavelength analysis, a gen- 
eral scenario for SNR 3C396 can be established as this: 
the remnant collides with a molecular wall in the west, 
flanking a southwestern pillar of molecular gas, and ex- 
pands in a low density region in the east (but encounter- 
ing a clump of backside molecular gas). 

The global dynamical evolution of the remnant can 
be inferred from the X-ray properties. Adopting kTx ~ 
0.9 keV as the average postshock temperature in the in- 
terclump medium, we have the mean blast shock veloc- 
ity v s = [16/cT x /(3/iraH)] 1/2 ~ 8.7 x 10 2 km s~\ where 
ttiyl is the hydrogen atom mass and ft = 0.61 is the av- 
erage atomic weight. Assuming an adiabatic expansion 
and adopting the SNR radius r s ~ 3'. 9 ~ 7.0^6.2 pc 
and the mean density of the interclump medium no ~ 

n x /4 - l/" 1/2 ^6.2 

cm , the SNR's age is estimated 
as t = (2r s )/(5v s ) ~ 3^6.2 kyr and the explosion en- 
ergy is E = (25/4^)(1.4n mH)^ 2 rf ~ 6(n /lcm -3 ) x 
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10 50 d 6 _ 2 erg (where £ = 2.026). The dynamical age esti- 
mated here is consistent to the ionization ages (the time 
since passage of the shock wave) of the X-ray emitting 
gas obtained from the spectral fit (Table EJ), except the 
region "W" for which the ionization age seems relatively 
low. 

4.2.2. Pressure discrepancy 

The thermal pressure of the X-ray emitting gas is about 

p h ~ 1 x 10- 8 /" 1/2 ^2 /2 d Y nes cm " 2 (Table Ej), which 
is about two order of magnitudes lower than that of the 
shock ram pressure p c \ ~ 3 x 10 -6 dynes cm -2 , as de- 
rived from the H2 emission of the Spitzer IRS observation 
towards the southwestern region (Hewitt et al. 2009). 
Similar pressure discrepancies have also been found in 
the cases of IC 443 (Moorhouse et al. 1991) and 3C 391 
(Reach & Rho 1999). Such discrepancies can be ex- 
plained if the western molecular gas is clumpy and there 
is a radiative shell between the hot gas and the dense 
molecular clumps (Chevalier 1999). A radiative shock 
has been revealed by the Spitzer IRS slit observation (He- 
witt et al. 2009). On the assumption of a crude pressure 
balance between the hot gas and the radiative shell, we 
have (Chevalier 1999) 



where n c \ is the hydrogen atom density in the clump, 
no the interclump density, and n rs the density in the 
radiative shell. The gas density in the dense clump, 
n c \ rsj 1 x 10 4 cm -3 , was estimated for the southwest- 
ern region from the H2 lines by two groups of authors 
(Lee et al. 2009; Hewitt et al. 2009). Considering that 
the X-ray emission represents the shock interaction with 
the interclump gas, we may assume no is about 1/4 the 
density of the hot gas (nx, see Table [3j) and thus is of 
order 1 cm -3 for the southwestern region. Hence we ob- 
tain n rs ~ 4 x 10 2 cm -3 . This number is similar to the 
electronic density (~ 3 x 10 2 cm -3 ) inferred from ionic 
lines (Hewitt et al. 2009). 

4.2.3. Supernova ejecta and progenitor 

In the north and south, the ~ 84kms _1 CO emission 
is very weak, the SNR does not seem to encounter the 
molecular wall, and the X-ray emitting gas is found to be 
metal enriched (Table EJ). We assume that the observed 
pattern of metal enrichment can reflect the composition 
of the supernova ejecta and thus be used to constrain the 
mass of the progenitor star (following, e.g., Rakowski et 
al. 2001; Chen et al. 2006). As the presence of a pulsar 
implies a core-collapse nature of the SN, we plot the rela- 
tive metal abundances predicted from such SN models for 
stars with normal metallicities (Z = Z@) (Thielemann, 
Nomoto, & Hashimoto 1996, hereafter TNH96; Woosley 
& Weaver 1995, hereafter WW95) in Figure dU The en- 
richment appears compatible with the WW95 model for 
a progenitor mass of M = 15M or the TNH96 model 
for a progenitor mass of M = 13M . A 13-15M pro- 
genitor is a B1-B2 star (Snow 1982) and ends its life as 
a type-IIP supernova (Heger et al. 2003). 

5. CONCLUSIONS 



We have investigated the molecular gas environment of 
the semicircular young composite SNR 3C 396 in multi- 
wavelengths, principally with 12 CO ( J =1-0 and J =2—1) 
and 13 CO (J=l-0) observations in SRAO and PMOD, 
and performed a Chandra spatially resolved thermal X- 
ray spectroscopic study of the SNR. The main conclu- 
sions are summarized as follows. 

1. We confirmed a cavity-like structure of molecular 
gas at - 67-72 krns" 1 in the 3C 396 SNR region, 
which was noticed by Lee et al. (2009) based on the 
GRS 13 CO (J =1-0) data. It is most probably to be 
located in the near distance (~ 4.3 kpc) according 
to the HI SA measurement. 

2. We find a series of evidence suggesting that the 
molecular gas at ~ 84 km s _1 is in physical con- 
tact with SNR 3C 396, which can reasonably ex- 
plain the multiwavelength properties of the rem- 
nant. Around this LSR velocity, a molecular wall 
appears to perfectly confine the western boundary 
of the SNR, which is bright in X-ray, IR, and radio 
emission. The CO emission fades out from west to 
east, indicating that the eastern region is of low gas 
density, accounting for the radio "blow-out" mor- 
phology in the east of the remnant. A molecular 
pillar is revealed in the SW, with one end intruding 
inside the SNR border. The shock interaction with 
this "pillar tip" can explain the X-ray and radio en- 
hancement in the SW and some infrared filaments 
there. The SNR-MC interaction is also favored by 
the relatively elevated 12 CO J =2-1 /J =1-0 line ra- 
tios at 80-83 krns -1 in the southwestern "tip" and 
the molecular patch on the northwestern boundary. 
Redshifted broadening (86-90 km s _1 ) is detected 
in the 12 CO spectral line profiles of the eastern 
81kms _1 molecular patch and may be the kine- 
matic evidence for shock-MC interaction. 

3. The establishment of the association between 
3C 396 and the 84 km s" 1 MCs results in a de- 
termination of the kinematic distance at ~ 6.2 kpc 
to the remnant, which agrees with a location at the 
tangent point along the LOS. 

4. The diffuse thermal X-ray emitting gas is found to 
be metal enriched except in the southwestern X- 
ray enhancement, and the metal enrichment in the 
north and south implies that the supernova pro- 
genitor is a 13-15M B1-B2 star. 

5. The X-ray spectral analysis suggests an age of the 
remnant of ~ 3 kyr. There is a pressure discrepacy 
between the X-ray emitting gas and the IR emitting 
cloud shock in the west, which can be explained 
with the inhomogeneous molecular gas and radia- 
tive shells between the hot gas and the dense molec- 
ular clumps. 
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TABLE 1 

Parameters of the ~ 84kms _1 Molecular Gas 



Region a 


AT(H 2 )(10 21 cm- 2 ) b 


M(M ) b 


Northwestern clump 


4.4/4.1 


1.9 x 10 3 d 2 2 /1.8 x 10 3 d 2 2 


Pillar 


4.6/4.5 


3.7 x 10 3 d 2 " 2 /3.6 x 10 3 d| 2 


Pillar tip 


4.3/4.4 


6.9 x 10 2 d 2 " 2 /7.1 x 10 2 c& 2 


Western wall 


4.2/3.8 


1.3 x 10 4 dg 2 /1.2 x 10 4 cZ|' 2 



a Regions are denned in Figure [8] 

b See text for the two methods used for the column density derivation. 
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TABLE 2 

Chandra X-ray Spectral Fits with 90% Confidence Ranges 



Region 


N 


W 


SW 


S 


E 


N+S 


AT H (10 22 cm- 2 ) 
/cTx(keV) 
Si 
S 

Ca 

r(10 n s cm- 3 ) 
fn e n x V/d 2 6 2 (10 57 cm- 3 ) 
X?(dof) 


4 84+ - 30 

^•°^-0.13 
, OO + 0.20 
i - OO -0.19 

1 66 +0 - 28 

-, 7 o+0.29 
— 0.27 

1 Qi+1-16 
i - yi -0.86 

2 19+ 1 - 50 

z - iy -0.63 

1 26 +0 - 48 

- LlZU -0.39 

1.27(78) 


5 50+ 97 

°- OU -0.92 
7 i+0.25 
U - ' ± -0.15 

92+ 31 

u.»__ 23 

1 99+ - 77 

- L - yy -0.67 
l a 

42+ - 46 
u -^ z -o.io 

9 7 o+4.32 
z - '°-2.73 

1.21(28) 


4 42+ - 41 

^•^ -0.39 

66 +0 - 10 

-, 9 o+0.23 
i._o_ 24 

1 20 +0 - 32 

j_._u_q 28 
l a 

>4.06 

o oc + 1-55 
z - oo -0.88 

0.66(48) 


4 74+0-26 

' -0.33 

1 04+ - 11 

1 42 +o.36 

l.60±g;f? 
3.08ti;2 

2 69 +2 - 79 

z - uy -1.09 
, 9 o + 0.52 
±._o_ 38 

1.40(59) 


4 95+ 1 10 
0.68_°;g 

l a 

l a 

l a 
>0.84 
2.34+g;£ 
1.14(19) 


4 82 +0 - 21 
^t.o__ 20 

, 91 +0.13 

1 - Z1 -0.13 

-L.OD_o.21 
i fi o+0.22 
l.DC)_ 2i 

1 Q 8 +0-93 

2 34+ 1 - 29 

Z - <34 -0.74 

2 50+0-71 

Z - OU -0.52 

1.54(101) 



a Fixed to the solar abundance. 



TABLE 3 
Parameters of the Hot Gas 



Region N W SW S E 

Volume(10 56 fd% 2 cm 3 ) a ' b 10.1 4.4 1.8 3.5 14.5 

n x(/" 1/2 ^2 /2 cm" 3 ) 1.0 2.3 3.3 1.7 1.2 

Mass(/ 1 / 2 cig / 2 2 M©) 1.2 1.2 0.7 0.7 2.0 

Ionization Age(10 3 / 1 / 2 4 / 2 2 yr) 5.6 0.5 > 3.2 4.0 > 1.9 

Pressure(10- 9 /- 1 / 2 (i^2 /2 dynes cm" 2 ) 6.1 7.3 9.7 8.1 3.5 



a Here / is the filling factor of the X-ray emitting gas and cfe.2 is the distance to the remnant in units of 6.2 kpc (Section 14.10 . 

b In the estimate of the volumes, we assume the oblate spheroids for elliptical regions N, W, SW, S, and E, with half axes / .88 x 1'.26 x l / .26, 

0'.49 x lMl x lMl, 0'.63 x 0'.63 x 0' .63, 0'.59 x 0'.91 x 0'.91, and l'.Ol x 1'.41 x 1'.41, respectively. 
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RA (J2000) RA (J2000) RA (J2000) 

Fig. 1. — SRAO 12 CO (J=l-0) intensity channel maps between 64 and 72 km s _1 (smoothed to a resolution of 0'.2 by interpolation), 
overlaid with the 1.0-7.0 keV X-ray (cyan, thin) contours and the 1.4 GHz radio (red, thick) continuum emission contours at levels of 1.5, 
3.5, 5.5, 7.5, 9.5 mJy beam -1 . 
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RA (J2000) RA (J2000) RA (J2000) 

Fig. 2. — SRAO 12 CO (J=l-0) intensity channel maps between 80 and 88 km s _1 (smoothed to a resolution of 0'.2 by interpolation), 
overlaid with the 1.0-7.0 keV X-ray (cyan, thin) contours and the 1.4 GHz radio (red, thick) continuum emission contours at the same 
levels as in Figure [T] 
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15.0 19:04:00.0 03:45.0 
Right ascension 



Fig. 3. Left: multiband image of SNR 3C 396. The PMOD 12 CO intensity map (> 2.2 K km s _1 or 11a) in 85-88 km s _1 interval is 
in red, the Spitzer 24 jum mid-IR emission in green, and the Chandra 1.0-7.0 keV X-rays (adapt ively smoothed with the CI AO program 
csmooth to achieve a S/N ratio of 3) in blue. Each of the intensity maps is in logarithmic scale. The five levels of the VLA 1.4 GHz 
contours are linear between the 78.6% and 12.4% maximum brightness. Right: Chandra ACIS X-ray tri-color image of SNR 3C 396. The 
Si (1.65-2.1 keV) is in red, the S (2.25-2.6 keV) in green, and the hard continuum (3.0-7.0 keV) in blue. 
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Fig. 4. The four left panels: the PMOD 12 CO intensity maps in the intervals of 67-72, 82-84, 84-86, and 87-90 km s _1 , smoothed 
to 2' to illustrate the large-scale confusion in t his direction. The four right panels: the PMOD 12 CO maps in the same velocity intervals 
after unsharp masking (see text in Section |3.1|) . overlaid by VLA 1.4 GHz contours at the same levels as in FigureO 
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80-83 km/s 83-85 km/s 85-88 km/s 
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Fig. 5. Maps of 12 CO J= 2-1/ J= 1-0 intensity ratio around the SNR 3C 396 at velocity ranges 66.5-68, 68-70.5, 70.5-72, 80-83, 
83-85, and 85-88 km s _1 after convolving the SRAO 12 CO (J =2-1) data to the same beam size as the PMOD 12 CO (J =1-0) data. 
Contours show the VLA 1.4 GHz emission at the same levels as in Figure 03 




70 80 90 70 80 90 

Velocity (km/s) Velocity (km/s) 

Fig. 6. — The CO spectra in 61-96 km s _1 interval for the eastern [the "East" region in Figure [4j 2.25 arcmin 2 centered at 
(19 h 04 m 14 s ,05°27 / 30 // )] and the southwestern [the "SW" region in Figure [1 2.25 arcmin 2 centered at (19 h 03 m 58 s , 05°25 / 30 // )] regions 
of SNR 3C 396. The thick black, thick blue, and thin lines are for the PMOD 12 CO (J=l-0), SRAO 12 CO (J=2-l), and PMOD 13 CO 
(J =1-0) (multiplied by a factor of 4) emission, respectively. The profiles of 13 CO ( J =1-0) and 12 CO ( J =2—1) are shifted by —1 and —2 K, 
respectively. The 12 CO (J =2-1) spectra are not convolved from beam size 48" to the 12 CO (J = 1-0) beam size 54". 
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3C396: 1420 MHz continuum 




20 40 60 50 100 20 40 60 50 100 

Velocity^ km /si Velocity(km/s) 

Fig. 7. — Upper left panel: the VLA 1.4 GHz continuum image of SNR 3C 396. Upper right panel: the HI spectra extracted from box 1. 
The HI spectrum of the bright source is from the solid-line box area, and the background spectrum is from the dashed-line box excluding 
the solid box area. The dashed line indicates 5 a derivation for velocities —100 ^ Vlsr ^ —50 km s _1 and 110^ Vlsr ^ 140 km s -1 . 
Lower panels: the VGPS HI (thin) and GRS 13 CO (thick) spectra extracted from boxes 2 and 3 (see the upper left panel), respectively. 
The vertical solid and dashed lines indicate the position of 69 and 84 km s -1 , respectively. 




04:25.0 19:04:05.0 45.0: 



8 10 12 14 16 18 20 22 24 26 

Fig. 8. PMOD 12 CO intensity map in 82-90 km s 1 interval, labelled with the defined regions for the parameter derivation of the 
84 km s _1 molecular gas. Contours show the 1.4 GHz radio continuum emission towards the remnant at the same levels as in Figure [3] 



Molecular environment and X-rays of SNR 3C 396 
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Fig. 9. — Chandra X-ray raw image of SNR 3C 396, labelled with the solid (dotted) defined regions for the source (background) X-ray 
spectrum extraction. The point sources (except the central pulsar) detected with a wavelet source-detection algorithm have been removed. 
The contours are the VLA 1.4 GHz radio continuum emission at the same levels as in Figure [3] 
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Molecular environment and X-rays of SNR 3C 396 
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Fig. 11. — The best-fit abundances of the X-ray emitting gas in the north and south (the 
compared with the predictions from core-collapsed SN models (TNH96 and WW95). 
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